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Planar superconducting junctions with a large effective Josephson coupling constant and a pro- 
nounced interface pair breaking are shown to represent weak links with small critical currents and 
strongly anharmonic current-phase relations. The supercurrent near Tc is described taking into ac- 
count the interface pair breaking as well as the current depairing and the Josephson coupling-induced 
pair breaking of arbitrary strengths. A new analytical expression for the anharmonic supercurrent, 
which is in excellent agreement with the numerical data presented, is obtained. In junctions with a 
large effective Josephson coupling constant and a pronounced interface pair breaking, the current- 
induced depairing is substantially enhanced in the vicinity of the interface thus having a crucial 
influence on the current-phase relation despite a small depairing in the bulk. 

PACS numbers: 74.50.-|-r, 74.20.De 



The Josephson current is one of the remarkable mani- 
festations of quantum coherence on the macroscopic scale 
in condensed matter physics. The supercurrent depends 
on the phase difference of the order parameters across 
the junction interface. The study of the current-phase 
relation (CPR) in the junctions makes it possible to iden- 
tify physical processes, which form supercurrents under 
diverse conditions. It is also beneficial for junction ap- 
plications. The problem attracted much attention while 
studying both highly transparent junctions with strongly 
anharmonic CPRs and tunnel junctions, where the sec- 
ond harmonic of the supercurrent comes into play due 
to the suppression of the first onei"— . The latter takes 
place in the junctions involving unconventional supercon- 
ductors with special interface-to-crystal orientations and 
at O-TT transitions. 

One of the earlier theoretical results, clarifying a vari- 
ety of aspects of the problem, is the anharmonic CPR for 
the superconducting point contactsii"— . Due to the neg- 
ligibly small current-induced pair breaking at any trans- 
parency value, the theory of point contacts is simpli- 
fied. The depairing plays an important role in form- 
ing anharmonic CPRs in highly transparent planar junc- 
tions, unlike its negligible role in point contacts. Since 
the critical current jc of usual planar junctions becomes, 
with increasing transparency, comparable with the de- 
pairing current j^p in the bulk, the junctions do not rep- 
resent weak links. In other words, in the junctions, the 
current-induced depairing brings about a pronounced an- 
harmonicity only in the crossover from the tunnel Joseph- 
son current to the bulk superflow^ii^"— . 

There are at least three types of pair breaking pro- 
cesses taking place in charge transport in the supercon- 
ducting junctions: the pair breaking produced by the 
phase-dependent Josephson coupling, by the current and 
by the interface itself. These are the very same effects 
which can lead to a noticeably anharmonic CPR not only 
at low or intermediate temperatures, but also near Tc- 
Here, I show that planar junctions with a large effective 
Josephson coupling constant and a pronounced interface 
pair breaking can possess strongly anharmonic CPRs and 



small critical currents satisfying the condition jc jdp- 
An enhancement of the current-induced depairing near 
the interface will be identified. The anharmonic super- 
current near Tc will be obtained within the Ginzburg- 
Landau (GL) theory in the presence of all three types of 
the pair breaking processes of arbitrary strengths. Along 
with the numerical solution based on GL equations, a 
new analytical CPR will be derived and shown to be in 
excellent agreement with the numerical data in a wide 
range of parameters. For tunnel junctions the obtained 
results present a new description of higher harmonics of 
the supercurrent and extend the known expressions for 
the first and second harmonics to include the effects of 
interface pair breaking. 

The CPRs obtained earlier near Tc with the micro- 
scopic boundary conditions for standard dirty s-wave 
junctionsi^ii^, have been considered in literature solely 
as the particular properties of the specific systems^. The 
anharmonic CPR obtained in this paper, and infiuenced 
by the interface pair breaking, is of general form inher- 
ent in the GL theory, and is applicable to a variety of 
planar junctions including those containing d^2_j^2-wave 
superconductors and/or magnetic interlayers. 

The free energy functional for Josephson junctions near 
Tc results in the GL equations and the boundary condi- 
tions (BC) for theni^Sr— . Consider symmetric junctions 
with a spatially constant width, which is much less than 
the Josephson penetration length, and with a plane in- 
terlayer at a; = of zero length within the GL approach. 
Assume the usual form of the GL free energy, which ap- 
plies, for example, to s-wave and d^2_y2-wa,ve junctions. 
If the Josephson coupling gj\'if+ — is strong, not 

only this term but all the interface and bulk contributions 
to the free energy generally participate in the formation 
of CPRs as a consequence of the dependence of absolute 
values of the order parameters at the interface on the 
phase difference. This concerns, in particular, the gradi- 
ent bulk term XjV^'p and the interface contribution of 
the form c/(|*+|^ + |*-|^). 

Moving on to the order parameter /(a:)e'^^^' normal- 
ized to / = 1 in the bulk without supcrflow, one gets the 
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first integral of tfie GL equation in the presence of the 
supercurren1j2^ in the form of 

Here x = x/S,, ^ = S,{T) is the superconducting coher- 
ence length, j is the spatially constant normalized cur- 
rent density j = j/jdp = -{^V2,/2){dx/ dx)p and foo is 
the asymptotic value of / in the depth of the supercon- 
ducting leads. 

The BC introduce in the GL theory at least two char- 
acteristic lengths £ = K/gj, 6 = K/g. The effective 
dimensionless Josephson g^ — gj^{T)/K and interface 
9s = gS,{T)/K coupling constants, associated with these 
lengths, will be used below. For symmetric junctions 
with / continuous through the interface, the BC for / as 
well as the expression for the Josephson current via /g 
and the phase difference x — X~ ~ X+ ^-t the interface, 
are obtained from the BC for complex order parameters: 

(2) 

Here the effective phase-dependent extrapolation length 

b{x) = {S~^ + 2i~^ sin^(^)) controls the pair breaking 
produced by the phase difference and by the interface. 
Let's denote gtix) = iss + 2g£sin^ ^). 

Since the material parameters in the normal state gj 
and g are not assumed to depend here on T near Tc, 
one should have \gi\ ^ 1 and/or \gs\ ^ 1 quite close to 
Tc due to large values of £,{T). However, the coupling 
constants gj and g can themselves be very small and 
the temperature range with large ge and/or gs be too 
narrow, as it occurs in standard tunnel junctions. Due 
to a very small surface pair breaking in conventional s- 
wave junctions, one parameter gg is usually assumed to 
describe the interfaces in ^ rather than both gg and 
gs as is in the regular case. At x = 0, such symmetric 
junctions contain no pair breaking at all, and the BC ([2]) 
is reduced to {df{0)/dx) = 0. 



Since only higher order terms begining with oc have 
been neglected in its derivation, the CPR (|4|) turns out to 
describe the current behavior almost perfectly, if j < 0.7. 
For j > 0.7 it gives a good interpolation of the numeri- 
cal solution based on resulting in the deviations not 
exceeding 10%. 

As seen in (jS]) and the anharmonic Josephson cur- 
rent j depends, in general, on the two dimensionless ef- 



If gj and/or g were very small, one would need to in- 
troduce into ([5]) the terms of the next order of smallness, 
in particular, in powers of the order parameter. Such 
terms could be of importance and bring about additional 
phase dependence and material-dependent parameters to 
the problem. Here, only the simplest conditions will be 
assumed, when ^ applies to a wide range of values 
of gi and gg. This agrees with the microscopic model 
rcsultsi^iiSi'^i^^"— and, for instance, takes place within 
the GL approach for sufficiently large values of gg and gs, 
which is the particular focus of this paper. 

There is no need to solve differential equation ([1]) in or- 
der to find /o, and, consequently, to find j via ([2]). One 
puts 2; = in ([T|) and, using eliminates the current 
and the first derivative of the order parameter. This re- 
sults in a biquadratic relation between the self-consistent 
values of /q and The second relation between them 
follows from the current conservation and the asymptotic 
formulas in the bulk. The current-induced depairing in 
the bulk is conveniently described via the superfiuid ve- 
locity j = {3V^/2)vs{l-v^s), fl = l-vj^. Equating 
the asymptotic expression for the current to that in ([2]) 
with /o = (1 — ^'s)c« one obtains Vg = ctgi sin x- Con- 
sidering that both quantities /o and foo as well as the 
current itself are now expressed via the only variable a, 
the fourth-order polinomial equation for a follows from 
the biquadratic relation between /o and foo 

2gUx)a - (1 - - a(« + 2)5,2 sin^ x] - 0. (3) 

Eq. ([3]) is exact within the conventional GL approach 
with BC ©. In the particular case of standard s-wave 
junctions, gb{x) = '^9esm^{x/2). Then ^ is reduced to 
Eq.(8) of Refji^, if one corrects a misprint Fs — ^ F^ 
in (8) and identifies the parameter of the GL theory 
g^^ = i/S, with the model parameter Tb entering the 
microscopic BC for dirty s-wave superconductors. 

An analytical solution of the problem can be obtained 
assuming a small depairing in the bulk <C 1 that allows 
to use ~ 1 ~ (4/27) j2 and to disregard the smaller 
terms on the right-hand side of (P). Then one gets from 
([1]), ([2]) a biquadratic equation for /o, which results in 
the analytical solution for the CPR: 



(4) 

I 

fective coupling constants gi, gs and the phase difference 
X- According to the simple physical arguments as well 
as the microscopic resulta^^"— , a variation of tunneling 
parameters principally modifies gi, while the surface pair 
breaking mostly contributes to gs ■ This signifies that the 
junction transparency D enters the combination of mi- 
croscopic parameters representing gg. The last statement 
agrees with the microscopic results for s-wave junctions 
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FIG. 1. a) jc as a function of ge, taken for various gg: 
1. gs = 0, 2. gs = 0.4, 3. gs — 1, 4. — 4. b) jc as a function of 
gs, taken for various ge: 1. ge = 0.1, 2. ge = 1, 3. ge = 100. 

with nonmagnetic interface a^^'^^i^^" — , where the corre- 
sponding combination is sometimes identified as the ef- 
fective transparencj«22i2^. The microscopic estimations 
of the effective Josephson couphng constant gi directly 
follow from those results. In the s-wave tunnel junctions 
{D < 1), one gets ge ^ D^{T){1-^ + ^~^), where I is 
the mean free path. In dirty superconductors, the ratio 
£,{T)/l can easily reach 100 even at low temperatures. 
Hence, for small and moderate transparencies, the quan- 
tity ge ^ D£^{T)/l can vary from vanishingly small values 
in the tunneling limit to those well exceeding 100 near 
Tc- In highly transparent junctions ((1 — D) ^ 1) the 
parameter ge oc (1 — D)~^ can be arbitrary largo^. The 
quantity ge can also take on negative values, which cor- 
respond to TT-j unctions, as seen in (|4]). 

The range of variation of the interface coupling gs 
can likewise be quite wide. For s-wave superconductor- 
insulator interfaces, the Josephson coupling vanishes and 
the extrapolation length b is reduced to d. The mi- 
croscopic estimations of S in such cases show it to be 
very large usually resulting in a negligibly small contribu- 
tion to the BC, unlike the superconductor-normal metal 
interfaces^. The length S can vary widely for d-wave 
superconductor-insulator flat surfaces, where it substan- 
tially depends on surface-to-crystal orientations^^"— . Al- 
though in this case 6 is strongly influenced by the surface 
roughness, in particular, by faceting^. 

A regular situation is characterized by a local suppres- 
sion of the order parameter at the interface. For this 
condition to hold, the effective extrapolation length b{x) 
should be positive at any phase difference and, hence, 
95 ^ (gs + '^gi) > 0. A superconducting state occurs lo- 
cally near the interface above the bulk Tc under the op- 
posite condition b{x) < with x ensuring the free energy 
minimunt^^. Only the simplest conditions gs/ > will be 
analyzed in detail in this paper, although the main results 
obtained here apply to substantially more general cir- 
cumstances. Other conditions, including magnetic field 
effects and/or negative gs,£, will be studied elsewhere. 

Figs. [T]a, b show the critical current jc as a function 
of coupling constants ge and gs- Solid curves have been 
calculated based on . Dashed curves correspond to the 
analytical expression Only for a small interface pair 



breaking {gs < 1) and for > 1, the current jc becomes 
comparable with 1, i.e., with the deparing current in the 
bulk. Thus the condition 5^ ^ 1 is the hallmark of a 
strong Josephson coupling. Comparatively small devia- 
tions of dashed curves from the solid ones are discernible 
only when the current exceeds about 0.7. With increas- 
ing gs , the growing interface pair breaking suppresses the 
critical current. For gs ^ 4, the critical current remains 
quite small jc <C 1 at any ge-, which would normally oc- 
cur in conventional tunnel junctions with small effective 
transparencies. In other words, in the regime of strong 
interface pair breaking gs > 4, the junctions represent 
weak links at any ge, including ge ^■ 

Though is a combined result of all depairing ef- 
fects, the origin of its characteristic anharmonic features 
is traced back unambiguously. The whole of the phase- 
dependence in except for that contained in gbix)-i is 
generated by the current via f^o on the right hand side 
or by the last term on the left hand side in ^ . Such de- 
pendence would retain the CPR (|4|) unchanged under the 
transformation X ~^ ^ ~ X- The symmetry is destroyed 
by the phase dependence of g^(x)i which originates from 
the BC ([2]) and can become pronounced, if \gs\ < '2\ge\- 

Whereas the CPR (|4]) is derived by assuming small 
depairing effects in the bulk, the depairing can be of cru- 
cial importance in Q within its domain of applicabil- 
ity. This is the case in the presence of a pronounced 
interface pair breaking, where an enhancement of the 
current-induced depairing, unlike the bulk, occurs near 
interfaces of junctions with ge 1. In particular, the 
phase-dependent term in the denominator in which is 
directly induced by the depairing, plays a key role in the 
case 3> 1 in restricting the normalized current value. 
The bracketed expression in the denominator originates 
from the coefficient before /q in the biquadratic equation 
for /q. The relative depairing correction coming from the 
bulk is (8/27)j^ = 2gjfQsm^x and its smallness signi- 
fies 2gjsm^xfo ^1- As seen, the term 2g'jsiv? x iii 
the denominator is allowed to exceed the unit consider- 
ably, when the condition 2gJ svc? xfo ^ ^ holds at the 
expense of a strongly suppressed order parameter at the 
interface /g ^ 1. Numerical results corroborate that, 
a gs ^4, the condition is satisfied at any ge including 
ge^ i (see also Figs. [T]a, b). This validates keeping ^ 
without its expanding in powers of gj sin^ x ^'^d explains 
the quantitative applicability of Q to junctions with the 
pronounced interface pair breaking at arbitrary ge- 

A number of specific CPRs follow from Q under a 
variety of particular conditions. Consider here two basic 
examples. The tunneling limit shows up in (j4|) under 
the condition \ge\ ^ 1. Developing (|4]) as series in ge at 
any value of gs, one obtains numerous harmonics whose 
weight is determined by ge and gs rather than by the 
transparency itself. The first and the second order terms 
result in 
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FIG. 2. a) CPRs j(x) for 9e = 5 and various gg: l.gs — 

4, 2.34 = 2, 3.gs = 1, 4.gs = 0.5 S.g* = 0. b) CPRs for 
gs = 4 taken for various gi: l.ge = 0.1, 2. gi = 1, 3. ge ~ 

5, -i.ge = 10, 5.ge = 50. 



Here j^'^ = (3V3/4).g^ (^2 + g| - Ig^l) is the main 

contribution to the first harmonic ji = jd sin x that 
is apphcable at any gg. Under the condition \gs\ <C 1 
it is reduced to the well-known result for tunnel junc- 
tions jco = (3-\/3/2)(7£, which is only justified when 
disregarding the interface pair breaking. In the oppo- 
site limit gg ^ 1 the pair breaking strongly suppresses 

the current and ~ jco/C^gs) ^ JcO: as is also 

known^i22r— i^i^^i^. In particular, the original current 

jii = ici^Jdp cc {T,-T) for 1 « 1 and j^l' cx {T,-T)^ 
for g^ ^ 1 near Tc. The second order terms in g^ bring 
about the main contribution to the second harmonic 
j2 = jc2 sin2x as well as corrections to the first one. The 
relative weight of the second harmonic in diminishes 
with increasing g^. The sign of jd coincides with the 
sign of gi , while the sign of jc2 is determined by the sign 
of gs. For small pair breaking < 55 ^ 1 the second or- 
der term cx g^ is simplified to the following correction to 
the current — -\/2.7co3f(sinx — (l/2)sin2x), in agreement 
with the corresponding microscopic resultsi^i^ for dirty 
and pure s-wave junctions. Note that the phase depen- 
dence generated by the current depairing shows up in (|4]) 
beginning with the third order terms in gi. 

The second example reveals the strongly anharmonic 
features contained in Consider junctions with the 
strong interface pair breaking ^ 1. Then a compara- 
tively simple approximate expression follows from (|4]) 

~ _ iV^gi sinx . s 

4[g| + 4(35 +5^)5. sin^f]- ^> 



The corresponding critical current jc = '^V^ge/^lgsigs + 
25^) I ^ 1 is always small. The associated phase dif- 
ference is determined by the relation sin^c = \gs{gs + 
'^gi)\/[{gs + gtY + gf]- it varies widely: Xc is smaU 
~ (gs/ge), if gi ^ gs^ and approaches 7r/2 in the op- 
posite limit gs 3> ge- Strongly anharmonic CPRs show 
up in ^ under the conditions gj ^ gg ^ 1. Also one 
has jc oc (Tc — T)^. Thus, at finite g, the temperature 
dependence jc{T) is quadratic quite close to Tc, where 
gs ^ 1- With increasing Tc — T, a crossover to the linear 
dependence on the temperature takes place in the region 
Tc - T <t: Tc, for sufficiently small g. 

Some of the CPRs j{x) are shown in Figs. [2]a, b. Ex- 
cept for the first curve in Fig. [5]b, the strongly anhar- 
monic CPRs in junctions with large Josephson couplings 
arc displayed. As seen in Fig. [2]a, the heights of the 
anharmonic peaks diminish considerably and the peak 
positions change weakly, when the interface pair break- 
ing goes up. Although the anharmonicity can be well 
pronounced even in the presence of quite a large pair 
breaking. This concerns, in particular, the curve 1 in 
Fig.[2]a, which is identical to the curve 3 in Fig.[3]b shown 
there in a different scale. Eq. describes the CPRs al- 
most perfectly and the corresponding dashed curves can 
be distinguished from the exact solid ones only near the 
high peak of curve 5 in Fig. [2]a. All curves in Fig. [2]b 
are also well approximated by a simple formula ([5]) with 
deviations (not shown) approaching only several percent. 
However, in contrast to ([5]) does not apply to describ- 
ing upper three curves in Fig. [5]a. The CPR similar to 
^ was found earlier within the microscopic description 
of the dirty s-wave junctions with metallic interlayersi^. 
The strong pair breaking can take place in those junc- 
tions, if the interlayer conductivity considerably exceeds 
the normal conductivity of the superconducting metal. 

In conclusion, the paper reveals the qualitative fea- 
tures and develops the quantitative description of the an- 
harmonic Josephson current near Tc- The interface pair 
breaking as well as the current depairing and the Joseph- 
son coupling-induced pair breaking have been taken into 
account and shown to play an important part in forming 
the CPR. The results obtained, in particular, concern 
the junctions involving d-wave superconductors and/or 
magnetic inter layers. 
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